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ABSTRACT 

We investigate the variation of the ratio of the equivalent widths of the FeIIA2600 
line to the MgIIAA2796,2803 doublet as a function of redshift in a large sample of 
absorption lines drawn from the Johns Hopkins University - Sloan Digital Sky Survey 
Absorption Line Catalog. We find that despite large scatter, the observed ratio shows 
a trend where the equivalent width ratio TZ = WFeii/WMgii decreases monotonically 
with increasing redshift z over the range 0.55 < z < 1.90. Selecting the subset of 
absorbers where the signal-to-noise ratio of the Mgll equivalent width tUMgii is >3 
and modeling the equivalent width ratio distribution as a gaussian, we find that the 
mean of the gaussian distribution varies as 77. oc (—0.045 ± 0.005)z. We discuss various 
possible reasons for the trend. A monotonic trend in the Fe/Mg abundance ratio 
is predicted by a simple model where the abundances of Mg and Fe in the absorbing 
clouds are assumed to be the result of supernova ejecta and where the cosmic evolution 
in the SNIa and core-collapse supernova rates is related to the cosmic star-formation 
rate. If the trend in 77 reflects the evolution in the abundances, then it is consistent 
with the predictions of the simple model. 

Key words: intergalactic medium; quasars: absorption lines; galaxies: halos, abun¬ 
dances, evolution 


1 INTRODUCTION 

The absorption lines of heavy elements (e.g., C, Mg, O, 
Fe) observed in the spectra of distant quasi-stellar ob¬ 
jects (QSOs) provide one of the best ways of observ¬ 
ing and understanding the evolving chemical enrichment 
of the Universe. Several lines of evidence suggest that 
these absorption lines arise in the circumgalactic medium 
(CGM) surrou nding galaxies lying close to the lines of sight 
to QSOs (e.g.. iB ahcbl fc SDitzeJ ll969l: lYannv et ^ Il990l : 


iBergeron fc Bois^ 199ll : ISteidel 19951 : Zhu et al. 2014 ). The 

source of the heavy elements in the CGM is not clear, but 
it is thought to result from enrichment of the halo gas 
by large-scale galactic outflows driven by a combination 
of stellar winds, supernovae, active galactic nuclei (AGN) 

* E-mail: dey@noao.edu 


and by the accretion of gas tidally st ripped during mergers 
from star-forming galactic disks (e.g., Op penl 2 eimeiyfc Davd 


l2006l : lD^alla Vecchia fc Schavel20^ : Suresh et aLll2015h . The 

ubiquity of galactic winds associated with star-formation in 
galaxies over a wide range in redshift is now empirically 
well establis hed (e-g., IWeiner et ahll ^OOgl: lErb et al.l l20ld : 
iMartin et al.|[2012 : iKornei et al.ll2oi2 ; Rubin et al. 2014 ) by 
the detection of blue-shifted absorption lines of Mgll and 
Fell in galaxies. These outflows can carry enriched material 
out of galactic disks to di s tances of > 50 kp c into galaxy 
halos teordoloi et al.l 20Uj;_ Rubin et al.l 2014 ) and thus en¬ 
rich the GGM (e.g., Prochaska et al.l 201ll) . Outflows of en¬ 
riched material from galaxy disks have also been invoked 
to regulate the mass-metallicity relationship observed in 


galaxies (e.g.. iTremonti et al.l 


IPave et al.ll2011 : Torrev et ahl 


200j; iFinlator fc Davdl2008l : 


20141 ): star-formation, fuelled 


by accreted gas, can increase a galaxy’s metallicity, but out- 
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flows can remove enriched gas from the galaxy and, in lower 
mass syste ms, prevent further s tar-formation and enrich¬ 
ment fe.g.. iMannucci et al.ll201(Ji . 

In this picture the metal enrichment of the circumgalac- 
tic medium is the result of supernovae and a natural conse¬ 
quence is that the ratio of a elements to Fe should reflect 
(to first order) the relative rates of core collapse supernovae 
(CCSN; which result from the explosion of a massive star) 
and Type la supernovae (SNIa; which result from the explo¬ 
sion of an accreting white dwarf). The heavy element yields 
from both cl asses of supernovae have been modeled in var - 
ious studies (iNomoto et alJ 1 19971 : iThielemann et ^ l2003ll , 
which typically And that the a-elements are formed pre¬ 
dominantly in CCSN, whereas Fe production is dominated 
by SNIa. The observed ratio of a-element to Fe equivalent 
widths in the spectra of absorption line systems over a range 
of redshifts can potentially provide an important window 
into the relative variation of SN rates, the efficiency of out¬ 
flows in enriching the CGM, and the resultant abundance 
patterns of diffuse CGM and IGM material with cosmic 
time. 

In this paper, we use the Johns Hopkins University 
- Sloan Digital Sk y Survey (JHU -SDSS) Metal Absorp¬ 
tion Line Gatalog (izhu &: MenardI l2013l . hereafter ZM13) 
to investigate the variation of the Fell/Mgll ratio as a 
function of redshift. The catalog tabulates the equivalent 
widths for 35,752 absorption line systems measured in the 
spectra of QSOs observ ed by the Sloan Digital Sky Sur¬ 
vey DR7 data release (lAbazaiian et al] l2009l : lYork et al.l 
I 2 OO 0 I I. We focus on the measurements of the equivalent 

widths of the MgIIAA2796.35,2803.53 (^S 3) and 

2 ’ 2 

FeIIAA2586.00,2600.73 {a^Dg —>■ z^Dr 9) transitions (see, 
e.g., IProchaska et al.ll 201 ll : iHartigan et al.iri999h . 

We present the sample in §2 and our analyses in §3 
and show that there is evidence for a trend of decreasing 
VUFeii/IUMgii with redshift. In §4 we discuss the possible 
origins of this trend. We consider the caveats and whether 
the ratio is a viable proxy for understanding the abundance 
ratio of Fe to Mg. We suggest that the observed trend re¬ 
flects the possible cosmic evolution of the Fe/Mg ratio and 
is consistent with the predictions of a simple model where 
the abundance ratio variations are driven by the ratio of the 
rates of SNIa to CGSN. 

Throughout this paper, we adopt a cosmology with 
[Dm, Da, crs, h] = [0.3,0.7,0.8,0.71 . Magnitudes are presented 
in AB units l|Oke fc Gunn lfl983l 'l unless otherwise specified. 


2 SAMPLE DEFINITION 


In this work, we m ake use of the JHU-SD SS Metal Absorp¬ 
tion Line Catalog l|Zhu fc Menardll 2013l. : hereaft er ZM13) 
based on the SDSS DR7 release l Abazaiian et al]l2009l 'FI. 

This catalog of 35,752 Mgll absorption line systems is 
constructed from 84,533 quasar spectra from SDSS DR7 and 
covers the redshift range between about 0.4 and 2.2 (de¬ 
termined by the wavelength coverage of the SDSS spectro¬ 
graphs; York et al. 2000). The window function, complete¬ 
ness and purity of the catalog were determined and cross- 


^ http://www.pha.jhu.edu/ gz323/jhusdss 


checked using different approaches, including Monte Carlo 
simulation for the whole sample and redshift path determi¬ 
nation for individual quasars (see ZM13 for details). The 
statistical incidence rate of the Mg II doublet as a function 
of redshift and strength was also presented in ZM13 and 
compared with the literature. 

The systematic effects in the catalog are also well un¬ 
derstood. The primary systematic effects are: (1) the im¬ 
perfect calibration of the relative spectrophotometry, due to 
imperfect F-star SED models; (2) the absorption induced 
by ISM/CGM of the Milky Way at given observer-frame 
wavelengths, most importantly at 3934A and 3969A (due to 
Ca H), 5891A and 5897A (due to Na I), and various dif¬ 
fuse interstellar bands; (3) poor subtraction of the telluric 
emission lines, which were, on average, under-subtracted in 
the reduction pipeline; and (4) relatively low completeness 
at both the short and long wavelength ends due to the 
reduced throughput at the blue end of the SDSS spectro¬ 
graphs and the poor subtraction of the bright OH lines at 
the red end. Finally, comparisons of the ZM13 catalog to 
an independently constructed DR7 catalog find good agree¬ 
ment in the measured equivalent widths in the two catalogs, 
with the difference in WmiiII measurement s less than 0.02 A 
l|Sevffert et akllioi^ : iGauthier et al.ll2014l L 

The ZM13 catalog therefore represents a well- 
understood dataset that spans a large range in redshift; 
the known systematics have negligible effects on any of the 
statistics derived here. The spectrophotometric calibration 
errors are at « 10 “® level (about an order of magnitude 
smaller than the photon Poisson noise) and therefore have 
no effect on the statistics presented here. The absorption 
by Galactic Ca II can be strong along some lines of sight 
and does affect the detection efficiency of extragalactic Mg 
H absorbers below redshift 0.43. We therefore choose not to 
include the absorbers at z < 0.43 in our sample. In addition, 
we restrict the fits and statisitical analyses to z > 0.55 to 
ensure all detected FeIIA2600 lines lie in the region of higher 
completeness at A > 4000A. The terrestrial lines are known 
to be under-subtracted and the noise (i.e., flux errors) at 
corresponding wavelengths were adjusted in accordance in 
the SDSS reduction pipeline. The Mg II detection pipeline 
has been fully-tested on cases with poor sky subtraction and 
shown to perform at the same efficiency as al ong other sight¬ 
lines. Thes e tests are discu ssed in detail in IZhu fc MenardI 
l|2013li and lLan et al.l (|2014| L At the red end, we restrict our 
analyses to z < 1.90, which corresponds to Mgll detections 
at roughly A < 8 IOOA. 

We extracted all measurements of the Mgll 
AA2796.35,2803.53, MgIA2853, Fell AA2586.00,2600.73 

absorption lines. The catalog contains a total of 35,752 
absorption line systems, >97% of which are “strong” 
absorbers with IUMgii 2796 > 0.3A. Since a few of these are 
only detected marginally, we restricted our analyses to a 
subset of 35,049 systems with absorption line signal-to-noise 
ratios of SNRMgii = WMgii/crwMgii > 3, where lUMgii is 
the sum of the equivalent widths of the Mgll 2796.35 and 
2803.53 resonance lines and o-WMgu i® the measurement 
uncertainty; 99.5% of these have IUMgii 2796 > 0.3. No 
restriction was placed on the strength of the Fell line 
(hereafter Fell refers to the 2600.73 line). 31,727 of the 
systems lie at 0.55 < z < 1.90 and constitute our primary 
statistical sample. 
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Figure 1. The histograms show sample distributions of the equivalent width ratios in four different redshift bins. The vertical dotted 
lines at VFFell/WMgll = 0.15 and 0.30 are shown as fiducials. The gaussian fits discussed in the text are shown as dashed curves with 
the peak value shown by the vertical dashed line. The legend in each panel lists the number of absorbers in the redshift bin (N), the 
value of the gaussian fit, the mean (/r) and standard-deviation (cr) of the best-fit gaussian, and the median redshift (< z >) of the 
absorbers within the bin. 


Some data in the catalog required editing. 2,554 sources 
have values of VKpeii < 10“®A, 816 of which lie at z < 0.55 
(and therefore outside our main sample) and the rest of 
which are distributed in redshift similarly with the over¬ 
all sample of absorbers. 7 more sources have (Twpen — 5A 
(most with meaningless values of many 1000s) and 3 have 
crwpeii ~ 0. In all these systems, we considered Fell to 
be undetected and, for the purposes of our analyses, we 
set Wpeii = O.OI 5 A, the minimum detectable equivalent 
width in the catalog, and the uncertainty o-Wp^n ~ 

In all cases, this resulted in crwpeii/IFpeii > 1.65. Within 
0.55 < z < 1.90, the redshift distribution of the < Icr Wpeii 
is very similar to that of the > la IFpeii equivalent widths. 


3 RESULTS AND ANALYSIS 

For each absorber, we defined the ratio TZ = RFeii/WAigii. 
The distribution of TZ shows a broad distribtion which is 
roughly gaussian with a weak tail to high values of TZ. Of 
the 31,727 absorbers, 4,965 (15.6%) have TZ < an- The red¬ 
shift distribution of these sources is similar to the overall 
distribution of absorbers. 

Given the large scatter in the points and the non¬ 
detections of Fell, we considered multiple approaches to de¬ 
termining whether there is a trend of TZ with redshift and 
estimating its significance. 

3.1 Fits to the VFpeii/IFMgii Point Distribution 

A least-absolute deviation fit to the overall point distribu¬ 
tion using 1000 bootstrap samples (with replacement) re¬ 
sults in a median slope of ATZjAz = —0.046 ± 0.003. The 
95% confidence limits on the slope are [-0.040,-0.052]. 

Given the significant fraction of Fell non-detections, we 
also measured the linear regression using the expectation- 
maximization (EM) technique and Schmitt’s binned linear 
regression method, as implemented in the ASURV package 
llLavallev et al.l[i992l : llsobe et al.lll986h . All points with TZ < 
a-jz were treated as upper limits. The EM method resulted in 
a slope estimate of —0.0406 ± 0.0035. Schmitt’s method was 
used with a redshift binning of S{z) = 0.15 and TZ binning 


of <5(77.) = 0.05 and resulted in a slope estimate of —0.034. 
We investigated different bin sizes, but the results did not 
change significantly. 

Finally, we used t he Bayesian approach to linear re¬ 
gression developed by iKellvl (l2007f) . This was very GPU¬ 
intensive and slow for our sample size, but resulted in a 
median slope of —0.0304 with a 95% confidence range of 
[-0.0343,-0.0252]. 

3.2 Fits to the Shape of the Wpeii/IFMgii 
Distribution 

We next investigate the redshift evolution in the shape of the 
distribution of 77 = lUpeii/IkMgii. In this approach, non¬ 
detections were just placed at their “measured” or edited 
value; given that the non-detections follow the redshift dis¬ 
tribution of the detected sample, any widening of the overall 
distribution is likely to be similar in all redshift bins. 

Following ZM13, we divided the sample into 
12 redshift bins defined between the redshifts 
z=[0.43,0.55,0.70,0.85,1.00,l.15,...,2.05,2.30] and mea¬ 
sured the distribution of the ratio 77 = Wpeii/IUMgii in 
each bin. We modeled the distribution in each redshift bin 
as a gaussian (fit using a non-linear least squares method) 
and found the best fit value for the gaussian mean. We also 
computed the median value of the distribution within each 
bin. Fig. [T] shows the distributions and the fits in a few 
sample redshift bins. We then performed linear fits (using 
least absolute deviation and least squares techniques) to 
the median and mean values to measure the variation of 77 
with redshift. 

In order to understand the uncertainty that may be in¬ 
herent in this analysis due to the effects of a complex and 
censored sample, we carried out a bootstrap analysis by se¬ 
lecting samples using replacement from the data in each red¬ 
shift bin. For each bootstrap sample, we carried out the same 
analyses (i.e., measuring medians, fitting gaussians, and fit¬ 
ting the median and means as a linear function of redshift). 

The results for the primary sample (i.e., 31,727 
W-Mgu/awMgu ^ 3 absorbers) are shown in Figure [2] The 
distribution of points is shown as a greyscale with the me¬ 
dian and quartiles represented as thin black lines. The best 
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fit peak and median values in each redshift bin are shown as 
large hlled triangles and circles, respectively, along with the 
uncertainties derived in the bootstrap analyses. The peak 
and median values of the TZ ratio distribution trend to lower 
values with increasing redshift. The same behavior is re¬ 
flected in the variation of the quartiles with redshift. A lin¬ 
ear fit to the peak values of the gaussian distribution within 
the redshift range 0.55 < z < 1.9 (i.e., using redshift bins 
with >900 absorbers) results in 7?. oc (—0.049 ±0.006)2 with 
95% confidence limits of [—0.055, —0.036]; a fit to the median 
values results in a similar slope of (—0.045 ± 0.005). 

The observed trend is not due to the censor¬ 
ship of the data. For our primary sample we used 
the least censored subset. However, we also investigated 
subsamples selected with a range of signal-to-noise ra¬ 
tios SNRFeiiorMgii > [0,1, 2, 3, 5, 7,10] and found significant 
trends (with TZ decreasing with increasing redshift in all 
cases) and with slopes (derived for hts to the median val¬ 
ues) ranging from —0.030 to —0.056, with a median de¬ 
rived slope value of « —0.045. The steepest (shallowest) 
slope values are found for the SNRMgii > 7 and SNRpeii > 0 
(SNRMgii > 10 and SNRpeii > 10) sample containing 20,813 
(2,962) absorbers (see Table 1). 

While higher SNR thresholds tend to select higher 
equivalent width lines, the correspondence between the two 
quantities is complicated and the equivalent width distri¬ 
bution at a given signal-to-noise threshold varies in a com¬ 
plicated way with wavelength. Hence, in order to further 
investigate whether the trend is dominated (or damped) by 
stronger absorbers, we investigated subsamples with a range 
of WmsU and Wpeii equivalent width cutoffs, selecting sam¬ 
ples with equivalent width cuts > [0,0.5,1.0,1.5,2.0] in either 
line. The derived slopes do not vary in a monotonic way with 
the equivalent width cuts, and there is no evidence from the 
current dataset that a subset of stronger or weaker absorbers 
are dominating the trend signal (see also Fig. fright panel). 

The trend of TZ with redshift appears to be primarily 
driven by an evolution in the Fell equivalent width. Fit¬ 
ting the medians of the distributions of the Mgll and Fell 
equivalent widths as a function of redshift results in mea¬ 
sured slopes of —0.06T0.05 for Wivigii vs. 2 and —0.20T0.02 
for Hpeii vs. 2 . The non-evolution in the slope of WAigii 
observed in our data is consistent wit h the results from 
previous high-spectra l-resolution studies dMshar et al.ll2007l : 
IChurchill et al.f[2003l i. 

3.3 Are Selection Effects Responsible for the 
Trend? 

We have established that the absorption line systems 
listed in the JHU-SDSS catalog exhibit a trend where the 
HFeii/VFMgii decreases with redshift. Although the catalog 
is constructed from the spectra of QSOs observed uniformly 
with the SDSS spectrographs, it is important to consider 
whether the observed trend could be the result of observa¬ 
tional selection effects. 

As an initial test, we investigated whether the line ratios 
Mgll A2796/MgH A2803 and Fell A2586/FeH A2600 in our 
sample show significant variation with redshift; they do not, 
and fits (to SNR>3 samples) result in slopes of —0.01T0.13 
and —0.03T0.01 respectively (see also the left two pan¬ 
els of Fig. [Sj). In addition, the MgIA2853/MgII equivalent 
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Figure 2. Distribution of the ratio of the Fell A2600 to Mgll 
AA2796,2803 equivalent widths as a function of redshift. The 
greyscale show the distribution of all absor ption line systems from 
the JHU-SDSS Absorption Line Catalog llZhu &: M4nardll201ljl 
with > 3cr WMgii measurements. The thin black lines show the 
median and quartiles of the distribution. The solid blue circles 
show the median values in different redshift bins; the solid pink 
triangles show the peak of gaussian fits to the ratio distribution 
(see text). The dashed lines show linear fits to the two sets of 
solid points. The grey hashed region denotes the range of slopes 
resulting from the Monte Carlo simulations of the null hypothesis 
(i.e., intrinsic TZ constant with redshift) described in the text. A 
small fraction (1.4%) of the absorbers have ratios larger than the 
y-limit shown. We include them in our analyses, but they have 
no effect on the results. 


width ratio does not vary significantly with redshift (slope of 
0.02T0.01). Changing the subsamples by varying the signal- 
to-noise ratio cutoff does not change these results signifi¬ 
cantly, although in the case of the Fell A2586/FeH A2600 
ratio, higher SNR cuts result in a slightly more significant 
detection of a trend (e.g., for SNR>5, Fell A2586/FeH A2600 
vs. redshift has a slope of —0.045 ±0.010). While the wave¬ 
length separation between these various features is smaller 
than that between the Mgll and Fell features, there does not 
appear to be any systematic effect introduced by the fact 
that the lines are measured in different parts of the spec¬ 
trum. The division in spectral coverage between the blue 
and red spectrographs at A ~ 5800A does result in a lower 
detection efficiency of absorption lines at these wavelengths 
and results in a slightly higher median equivalent width for 
the 2 « 1.2 bin; however, there is no obvious feature in the 
distribution of the line ratios and this bin does not affect 
the determination of the overall trend. 

In order to determine whether the complicated selec¬ 
tion function for absorption line detection could be respon¬ 
sible for the trend, we used the measured selection function 
from ZM13 to create 1000 simulated datasets as follows. 
We simulated a large number of Mgll absorption lines (dis¬ 
tributed randomly in redshift within each redshift bin) using 
the best-fit intrinsic exponential distribution of rest equiva¬ 
lent width computed by ZM13 (see equation 3 and table 1 
in ZM13). For each Mgll absorber, the HMgii 2796 /IFMgii 2803 
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SNR(Mgll) 

>0 

>1 

>2 

SNR(Fell) 

>3 

>5 

>7 

>10 

>0 

-0.0430±0.0048 

-0.0468±0.0044 

32167 

-0.0506±0.0046 

-0.0678±0.0034 

27042 

-0.0590±0.0044 

-0.0569±0.0030 

20487 

-0.0522±0.0060 

-0.0522±0.0027 

14873 

-0.0417±0.0046 

-0.0432±0.0030 

8687 

-0.0355±0.0061 

-0.0372±0.0033 

5483 

-0.0325±0.0069 

-0.0339±0.0043 

3044 

>1 

-0.0428±0.0050 

-0.0470±0.0043 

32157 

-0.0505±0.0047 

-0.0674±0.0034 

27034 

-0.0586±0.0044 

-0.0570±0.0030 

20483 

-0.0524±0.0056 

-0.0525±0.0027 

14869 

-0.0417±0.0046 
-0.0431±0.0030 

8686 

-0.0353±0.0062 

-0.0373±0.0034 

5483 

-0.0331±0.0071 

-0.0339±0.0042 

3044 

>2 

-0.0432±0.0048 

-0.0468±0.0042 

32137 

-0.0508±0.0047 

-0.0673±0.0035 

27021 

-0.0590±0.0045 

-0.0569±0.0028 

20474 

-0.0527±0.0056 

-0.0525±0.0028 

14866 

-0.0420±0.0047 

-0.0434±0.0031 

8686 

-0.0351±0.0063 

-0.0372±0.0034 

5483 

-0.0325±0.0066 

-0.0338±0.0041 

3044 

>3 

-0.0445±0.0049 

-0.0487±0.0045 

31727 

-0.0496±0.0046 

-0.0667±0.0034 

26864 

-0.0583±0.0042 

-0.0568±0.0030 

20447 

-0.0516±0.0053 

-0.0524±0.0028 

14851 

-0.0417±0.0047 

-0.0432±0.0030 

8679 

-0.0351±0.0061 

-0.0370±0.0034 

5478 

-0.0319±0.0068 

-0.0336±0.0040 

3043 

>5 

-0.0543±0.0051 

-0.0582±0.0042 

27602 

-0.0478±0.0045 

-0.0634±0.0036 

24397 

-0.0521±0.0047 

-0.0525±0.0028 

19655 

-0.0503±0.0056 

-0.0520±0.0028 

14663 

-0.0410±0.0046 

-0.0430±0.0030 

8647 

-0.0333±0.0061 

-0.0370±0.0034 

5457 

-0.0315±0.0069 

-0.0336±0.0043 

3035 

>7 

-0.0578±0.0043 

-0.0608±0.0037 

20813 

-0.0487±0.0054 

-0.0599±0.0034 

19357 

-0.0431±0.0045 

-0.0473±0.0030 

16896 

-0.0481±0.0048 

-0.0474±0.0026 

13619 

-0.0416±0.0045 

-0.0438±0.0029 

8495 

-0.0337±0.0059 

-0.0372±0.0035 

5396 

-0.0308±0.0069 

-0.0340±0.0042 

3008 

>10 

-0.0563±0.0078 

-0.0610±0.0039 

13483 

-0.0499±0.0064 

-0.0601±0.0039 

13028 

-0.0459±0.0048 

-0.0491±0.0033 

12234 

-0.0441±0.0042 

-0.0455±0.0029 

10879 

-0.0427±0.0051 

-0.0448±0.0030 

7859 

-0.0364±0.0054 

-0.0410±0.0032 

5222 

-0.0305±0.0067 

-0.0353±0.0043 

2962 


Table 1. Trend slopes for subsamples defined by censoring the data by the signal-to-noise ratio in Wpeii and ITiyigii- For each pair 
of signal-to-noise ratio cuts, the top entry is the median slope and Itr error (determined using bootstrap resampling) determined by 
fitting (using least absolute deviation) the median values In each redshift bin; the middle entry is the median and la error for the 
slope determined by a linear least squares fit to the mean of the gaussian model for the distribution; and the bottom entry is the total 
number of absorbers in the censored sample. The primary sample used In this paper is the one listed in bold face for SNR(Mg)>3 and 
SNR(Fe)>0. 


doublet ratio was drawn from a lognormal model based on 
a fit to the observed distribution of the ratio in the JHU- 
SDSS dataset. Then, for each Mgll absorber, an Fell ab¬ 
sorption line was drawn randomly from a normal distribu¬ 
tion of TZ with a mean of 0.233 and cr=0.143. In order to 
test whether the selection would generate a systematic bias, 
the TZ distribution was assumed to be constant with red¬ 
shift. Equivalent width uncertainties for the Mgll and Fell 
lines were assigned randomly, drawn from model fits to the 
actual fractional error distributions in these two lines. In 
each case, we constructed the two-dimensional distribution 
of aw/W versus W and modeled the surface as a lognormal 
in aw/W and a polynomial in W. We then “observed” this 
set of intrinsic Mgll+Fell absorbers using the equivalent 
width selection function for the JHU-SDSS catalog (ZM13). 
We selected absorbers and fit each of the 1000 such datasets 
in the same manner in which we fit the observed catalog and 
derived slope measurements for any gradient in WFeii/WMgii 
ratio with redshift. Since there was no evolution assumed in 
the data, any resulting slope is purely due to selection ef¬ 
fects (i.e., completeness as a function of equivalent width 
and observed wavelength) and the variation in the intrinsic 
distribution of WMgii with redshift. 

We do not detect the observed slope in the simulated 
TZ distribution with a high degree of reliability. The 1000 


trials resulted in a normal distribution of slopes with a small 
negative bias (median recovered slope = —0.0148; standard 
deviation of 0.0033); the 95% (99%) confidence limits on the 
slope are [—0.021, —0.008] ([—0.023, —0.0062]) and none of 
the 1000 trials resulted in slope determinations < —0.0245. 
We therefore conclude that the observed trend is not the 
result of selection effects. 


4 DISCUSSION 

Our results show that the ratio of the equivalent widths 
of Fell to Mgll absorption lines along the lines of sight to 
QSOs decreases with redshift. The trend survives a range of 
signal-to-noise ratio or equivalent width cuts and does not 
appear to be due to selection effects within the JHU-SDSS 
ZM13 dataset. In this section, we discuss possible physical 
origins for the equivalent width ratio evolution. 

The strength of an isolated absorption line as mea¬ 
sured by its equivalent width is related to its ab undance 
in th e gas phase and the veloc ity spread of the gas (iMenzell 
1 19361 : iBakeil [l936l : I Allerj [l943 ). It is tempting to interpret 
the WFeii/IUMgii ratio in terms of the Fe/Mg abundance ra¬ 
tio. However, in order to interpret the observed equivalent 
widths in terms of physical quantities associated with the 
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absorbing clouds, we need to consider the effects of (i) sat¬ 
uration, (ii) kinematics, (iii) ionization, and (iv) depletion. 
We consider each of these in turn. 

(i) Line Saturation: The WMgii distribution for the SNR>3 
sample is a lognormal which peaks at « 2A, implying that 
the individual doublet lines have typical equivalent widths 
>lA and lie in the logarithmic (i.e., “saturated”) part of 
the curve of growth. It is therefore not obvious that we can 
interpret any observed trend with redshift as resulting from 
a trend in the abundance ratio of Mg to Fe without first 
understanding the effects of line saturation. 

In order to investigate the effects of saturation on the 
observed ratio, we computed the curves of growth for Mgll 
and Fell for two different values of the Doppler parameter, 
6 _d = 10 km s“^ and 50 km s“^. For each case, we computed 
the expected TZ values that result from a wide range in Mgll 
column density 12 < log(A^Mgii) < 17.9 and [Mg/Fe] varying 
between 1 and 10 times the solar value. For a given bo, even 
strongly saturated lines result in the equivalent width ratio 
TZ depending on the [Mg/Fe] abundance ratio. Hence, even 
if the observed absorption lines are single and saturated, 
they only dampen, but not erase, an intrinsic variation of 
the abundance ratio. Most Mgll absorptio n systems with 
Wmb u > O.SA are multiple systems (e.g., IChurchill et al.l 
120030 . and this complicates the simple picture described 
here. Fully modeling the result of saturation and multiplicity 
is beyond the scope of this paper, but instead we look to the 
data for clues to whether saturation effects are responsible 
for the observed trends. 

Fig. [3] shows the effect of saturation on three pairs of 
equivalent width ratios in our observational dataset. The 
IFMgii2796/IFMgii2803 and WFeii2586/H4’eii2600 ratios (sliown 
in the first two panels) asymptote toward the saturation 
value with increasing VKiugii. However, neither of these ra¬ 
tios shows any evolution with redshift. In contrast, the 
ITFeii/VFMgii ratio clearly evolves, with the median line ratio 
decreasing with increasing redshift in all VTiugii bins. Fig. [3] 
also shows that the higher equivalent width bins also exhibit 
the evolution with redshift, albeit at lower strength than 
the lower equivalent width (and less saturated) bins. This is 
consistent with our expectation that saturation effects will 
tend to weaken the relationship between TZ and the intrinsic 
Fell/Mgll abundance ratio of the absorbers, and then that 
these effects are likely only to dampen (rather than induce) 
a potential trend in TZ with redshift. 

(ii) Multiplicity & Kinematics: 

If the observed absorption lines are instead the result of 
blends of multiple unresolved (lower equivalent width) lines, 
then the individual lines might arise in lower column density 
clouds and are therefore presumably less saturated. The real 
situation is likely to be more complicated, with lines forming 
in clouds with varying kinematic properties. 

High-resolution studies have demonstrated that Mgll 
lines are often made up of a number of components, typically 
between 2 and 1 0, with the number of components increas¬ 
ing with HAigii ilChurchill et al.ll2003l l. However, they also 
find that the Mgll and Fell absorption lines show similar 
Doppler parameters and correlated column densities and the 
distribution of A^(FeII)/A'^(MgH) column ratios does no t de¬ 
pend strongly on velocity spread dChurchill et al.ll2003l l , nor 
does the kinematic spre ad show any eviden c e for significant 
evolution with redshift dMshar et al.l l2007l : IChurchill et al.l 


l200(tl . Hence, assuming the Fell and Mgll absorption lines 
are equally sensitive to the variations in the clouds as they 
evolve, the evolution of the ratio WFeii/HMgii for individ¬ 
ual absorption systems should change monotonically with 
the abundance ratio A^eii/A^Mgii, although the magnitude 
of the change may be dampened by saturation effects. If, for 
some reason, the velocity structure or Doppler parameter 
for these two species differs, due to e.g., differences in ion¬ 
ization structure or incomplete mixing of metal species, this 
monotonic relation can break down. For example, incom¬ 
plete metal mixing could arise in a delayed Fe-enrichment 
scenario: if the kinematic spread of Fe-rich material is nar¬ 
rower than that of a-enriched gas, this could give rise to 
a suppressed WFeii/WMgii ratio without it reflecting a true 
dearth of Fe atoms. However, even if this somewhat con¬ 
trived scenario were realized, it would still point to an evo¬ 
lution in the sources of CGM enrichment. While it may be 
possible to create a trend in TZ with redshift as a result of 
an evolving bn, there is no obvious reason why the Doppler 
parameter would evolve with redshift for just one species. 
Better spectroscopic data (higher spectral resolution, and 
on more species and ionization states) on much larger sam¬ 
ples than currently exist can help explore this issue further. 

(iii) Ionization Correction: The first two ionization po¬ 
tentials for Mg and Fe are very similar: for Mgl^Mgll 
and Fel—>-FeH, they are 7.646 eV and 7.902eV respec¬ 
tively; the second ionization potentials (i.e., Mgll^Mglll 
and Fell—>FeI II) ar e 15.0357 eV and 16.1878 eV respec¬ 
tively llMortorJl2003l : iKramida et al.ll2013ll . The ratio of the 
ionic abundances of the two elements, A^Mgii/Aipeii, is de¬ 
termined by the shape (at these energies) of the spectrum 
of the ionizing metagalactic radiation, which does not vary 
sig nificantly over the redsh ift range considered here. Using 
the lHaardt fc Madaul (l2012f) model for the UV metagalactic 
flux as a function of redshift, we find that the ratio of the 
emissivity at the Mg and Fe 1st ionization potentials varies 
by less than 4% for 0.4 < z < 2.0; the emissivity ratio at 
the wavelengths of the 2nd ionization potential varies by less 
than 2%. 

The bulk of the absorbers (>96%) in the JHU-SDSS 
Metal Absorption Line Catalog have IFMgii2796 > 0.5A. 
These absorbers therefore correspond to regions with hy¬ 
drogen column densities of log(A Hi) > 19.4 (if the 

Dopp l er parameter is 10 k m s“^; iMenard fc Chelouchd 
120091) . ICiavalisco et al.l ll201ll ') calculate the variation of the 
N (Mgll) /N(Fell) ratio with ionization parameter for the 
iHaardt fc Madaul lll996ll metagalactic radiation field. They 
find that for hydrogen column densities log(A^Hi/cm“^) > 
19, the N(MgII)/N(FeII) ratio (and hence its inverse) is es¬ 
sentially constant over a wide range of ionization parameter 
(see also ICrighton et ^l2015h . 

It is therefore unlikely that differences in ionization 
could drive the observed trend. 

(iv) Depletion: Mgll absorption line systems are known 
to contain dust, exhibiting dust-to-ga s ratios comparable 
to th at observed for normal galaxies llMenard fc Fukugital 
I 2 OI 2 II . Since both Mg and Fe can be depleted onto dust 
grains, it is necessary to understand whether differing de¬ 
pletion rates for the two elements as a function of redshift 
can result in the trend observed in their equivalent width 
ratios. 

Based on studies of the interstellar medium in our 
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Figure 3. The impact of line saturation, shown as the variation of VFMgll 2796 /W ^2803 (Left panel), WFell2586/W^Fe2600 (Middle panel), 
and Wpeii/W^ivigll (Right panel) with the Mgll equivalent width. The small points show the distributions for the > 3 (t sample described 
in the text and the large coloured circles show the median values of the distributions in the redshift bins shown in the legend. The left 
two panels show that the ratio asymptotes to a constant value with increasing Wniglli suggesting that the lines are becoming saturated, 
but that the behavior is the same at all redshifts. In contrast, the Wpeii/R^Mgll (Right panel) shows a systematic evolution with redshift 
in all IFMgll bins. 


Galaxy, the depletion rate for Fe is approximately 2—7 times 
larger than that for Mg (e.g.. ISavaee &: Sembachlll996l ). i.e., 
Fe is depleted from the gas phase at a higher rate than is 
Mg. Using the tabulated model of Ijenkinj ll2009l) for the 
dependence of the Mg and Fe gas phase abundance on the 
line-of-sight depletion strength F ,, we find that the column 
density ratio log(A^Mg/A^Fe) varies as (0.29±0.06)F’*, i.e., the 
gas-phase column density ratio increases with increasing de¬ 
pletion and dustier absorbing clouds (which would presum¬ 
ably be more depleted in Mg and Fe) would exhibit higher 
values of log(A^Mg/A^Fe)- Hence, the only way to reproduce 
a trend where the WFeii/LFMgii ratio decreases with redshift 
through purely dust depletion effects is to have higher red¬ 
shift absorbers be dustier (and more depleted) than lower 
redshift absorbers, an unlikely scenario. 

In contrast, we expect Fe and dust content to increase as 
the Universe ages and that Fe is likely to be fractionally more 
depleted onto dust at lower redshift. Hence, any detection 
of an increase (to lower redshift) in the observed gas-phase 
Fe/Mg abundance ratio is likely to reflect a stronger intrinsic 
increase in the abundance ratio. In other words, a measure¬ 
ment of the gas-phase Fe/Mg ratio would underestimate the 
true Fe/Mg, and the factor by which it is underestimated 
would be larger at low redshift than at high redshift. This 
goes in the sense, therefore, of making the observed relation 
imply a stronger underlying variation. It is therefore unlikely 
that depletion effects could explain away the observed trend. 


4.1 A Simple Physical Model 

In contrast to various possible (and somewhat contrived) 
explanations discussed above, the expected variation of the 
rates of SNIa and CCSN should naturally produce the ob¬ 
served variation in TZ and we therefore consider this next. 

We model the average universal Fe/Mg abundance ra¬ 
tio as resulting from the supernovae formed at each epoch 


as a result of the cosmic star formation rate history^ We 
assume that the rate of core collapse supernovae follows 
the c osmic star formation rate density, i.e., rccsN oc ipiz) 
(e.g.. [Madau fc DickinsonI l2014h . Since CCSN progenitors 
are massive (M > 8 Mq) stars, the rate of CCSN is closely 
tied to the cosmic star-formation rate density in the Uni¬ 
verse, i.e., 

f?cc(.z) = 0.0068^^(2) yr~^ Mpc~^, ( 1 ) 


where is the cosmic s tar formation rate densit y 

(in units of MQyr“^Mpc“®; iMadau fc DickinsorJ [20l3 ). 
Rcc{z) assu mes a Salpeter initial mass function with 8 < 
M < 40 Mm llMadau fc DickinsorJ l2014h . In contrast, SNIa 
rates exhibit a “delay time distributio n” (DTD) relative 
to the time since s tar fo r mation (e.g ., Totmii^^_alJ_200^ 
Mjaz Sz Mannu^ I 2 OI 2 I : iMaoz et al.l I 2 OI 2 I : IClaevs et al.l 


2014h . Since there is much debate about the exact form of 


this distribut ion, here we assume the si mple power-law form 
suggested bv iMaoz fc Mannuccil (l2012l ): 


Ria.{t) = f '0(1 — 7")DTD(r)dryr ^ Mpc (2) 

Jo 

where DTD(r) = 4 x 10~^®(r/lGyr) yr“^Mg^. 

We compute the evolving mass density of each element 
as the sum of the SN yields: 

Mx~Yx^^^ f Rsm3.it)dt+ Yx‘^^^ f i?ccsN(t)dt, (3) 
Jo Jo 

where Tx is the average mass yield for metal X (in Mq) for a 
SN of type S. The number density of element X is MxjAx, 
where Ax is the relative atomic mass of element X. The 
SN yiel ds are very uncertain: based on a review of the lit¬ 
erature, IWiersma et al.l (l2009l ) conclude that calculated nu- 
cleosynthetic yields from SNe vary by factors of ~2. Here, 


^ We ignore the contribution of AGB stars here because (a) su¬ 
pernovae dominate the metal enrichment, and (b) we are con¬ 
sidering absorption line systems which likely are clouds at large 
distances from the galaxy nucleus or disk, where the enrichment 
is probably due to energetic SN ejecta. 
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Redshift 

Figure 4. The blue dashed line shows the variation of Fe/Mg 
abundance ratio with redshift predicted by the simple model 
described in the text, based on a CCSN rate tied to the 
Madau &: Dickinson (201 4) star-formation rate density and the 
iMaoz fc Mannuccil ll2012l~l power-law delay time distribution for 
the SNIa rate. The dot-dashed purple line s hows the variation 
with redshift of the Fe/Mg ratio computed bv lYates et al.l <l2013l') 
for a model Milky Way - type galaxy. The red circle with the 
error bar represents the solar abundance ratio. The solid line and 
the hashed region surrounding it show, respectively, the observed 
variation in the equivalent width ratio TZ and the 99.5% confi¬ 
dence limits on the slope. We have normalized the fit to match 
the solar value at 2 = 0, under the assumption that the abun¬ 
dances are linearly proportional to the equivalent widths. 


we assume that the typical SNIa produces 8.6 x 
of Mg and O.74M0 of Fe, and that the typi cal CCSN pro- 
duces 0.12 Mq of Mg and 0.09 M© of Fe l|Nomoto et al.l 
[1993). Choosing a single number for the average yields is a 
simplification, but our aim here is to explore whether the 
overall variation observed in the VfFeii/WMgii ratio can be 
explained by just the simple evolution of the SN rates. 

The results are shown in Fig. [I] The solar abundance 
ratio of (Vpe/VMg)© = 0.794lg'Jg:!. is also shown (based on 
the Mg and Fe photosp heric abundances of 7.60T0.04 and 
7.50±0.04 respectively: lAsplund et ^ l2009li . This simple 
model over-predicts the observed SNIa rates by a factor of 
2 and over-predicts the solar photospheric abundance ratio 
of Fe/Mg at redshift a = 0 by « 40%. This is not surprising 
given the large uncertainties in the yields and our other as¬ 
sumptions (e.g., ignoring the contribution from AGB stars, 
uncertainties/variations in the IMF, DTD, SFRD, and the 
large scatter i n the data)n 

Recently, lYates et al.l (l2013l ) modeled the metal enrich¬ 
ment in a galaxy resulting from AGB stars and SNe, in¬ 
vestigating different assumptions for SN yields and SNIa 
DTDs. They find that the observed local mass-metallicity 
relationships and the [0/Fe] abundance in the Milky Way 
can be reproduced by models in which the SNIa DTD is a 
power-law and the light a elements are expelled in galactic 
wi nds. The redshift v ariation of the Fe/Mg ratio computed 
bv lYates et al.l (l2013l f for a model Milky Way - type galaxy 
using a power-law DTD is also shown in Fig. 2] rather re- 


® Correcting the RsNia rate by a factor of 0.6 can result in both 
a better fit to the observed rates and the solar abundance ratio. 


markably, its slope matches the evolution in TZ observed in 
the JHU-SDSS dataset. This lends further, though still cir¬ 
cumstantial, support that the observed evolution can be ex¬ 
plained as the result of the increasing dominance of SNIa 
yields changing the mean composition of the outflowing gas 
responsible for the absorption line systems over time. 

The comparisons of the models to the slopes measured 
in the data are only meant to be illustrative and not dis¬ 
criminatory. The scatter in the slope measurement is large 
and the models are only notional at this stage. Nevertheless 
it is instructive that the measured variation shows slopes 
comparable to those predicted by models where all the evo¬ 
lution in the abundance ratios is the result of star formation 
in galaxies and enriched material ejected into their gaseos 
halos. The large scatter in the data may reflect the diversity 
of star formation histories, impact parameters of QSO lines 
of sight through halos, and inhomogeneities in the halo gas. 

Individual absorbers at any given redshift show signif¬ 
icant scatter and the trend (of TZ with z) reported here is 
subtle; it was only possible to measure the trend using the 
very large catalog of absorbers resulting from the SDSS QSO 
spectra. It is therefore not surprising that previous high- 
resolution studies (e.g ., ChurchilLet_uL 200^Msharet_3 


I 2 OO 7 I : [Narayanan et aluOOSl : Rodriguez Hidalgo et al.ll2012l l 

did not see the trend of TZ with redshift in the s mall sam¬ 
ples (< 100 absorbers) st udied (however, please seelSrianandl 
(Il99dl ’s analysis of the ISteidel &: Sargent ll 19921 ) dataset, 
where a possible trend of decreasing WpeiiA 2382 /kFiUgii is re¬ 
ported) . In a VLT/UVES high-resolution study of 100 weak 
(i.e., IFMgii 2796 < O. 3 A) Mgll absorbers, iNaravanan et al.l 
( 20081) found that clouds with high ITFeii/VFMgii disappear 
from their sample at high redshift {z > 1.2). In a com¬ 
plementary st udy of 81 strong absorbe r s iden tified in the 
same dataset, iRodrfguez Hidalgo et al] (|2012|) found that 
the WFeii/VFMgii ratio for IFMgii2796 > lA absorbers is typi¬ 
cally high and exhibits some evolution with redshift; weaker 
absorbers (0.3A< IFMgii 2796 < lA) instead showed a wider 
range and scatter in the IFFeii/WMgii ratio with no obvious 
trend. The authors interpret this as primarily due to kine¬ 
matics: high-redshift absorbers with IFMgii 2796 > lA tend 
to be formed in regions with more kinematic spread (and 
the Fell and Mgll lines are less saturated) than their low- 
redshift counterparts. They also consider the effects of ion¬ 
ization and a-enhancement, but conclude that they cannot 
discriminate between these effects. 

The SDSS spectra have insufficient resolu tion to ex¬ 
plore this issue adequately (lYork et al.l I 2 OOOI ). However, 
splitting the sample into two equivalent width bins — 
0.3A< IFMgii 2796 < lA and WMg ii2796 > lA fol¬ 
lowing IRodrfguez Hidalgo et al] (l2012f) — we find that 
both subsets show evolution (with slopes and 95% confi¬ 
dence ranges of -0.0502, [-0.0591, -0.0422] and -0.0434, 
[—0.0516, —0.0344] respectively). There is also evidence that 
the TZ values show a much broader distribution for the 
0.3A< IFMgii 2796 < lA subset than for the IFMgii 2796 > lA 
subset. It is not obvious from the large SDSS dataset, there¬ 
fore, whether splitting the data by equivalent width is jus¬ 
tified or if the absorber population should be treated as a 
continuum, reflecting a range of gas physical conditions in 
the halos of galaxies. 

Given the various caveats already described in using 
low-resolution data, it will be necessary to obtain high reso- 
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lution spectroscopic observations of a sufficiently iarge sam¬ 
ple 1000s) of absorbers to properly disentangle the effects 
of saturation, kinematics and abundance. A more detailed 
investigation of the scatter in the equivalent width and ratio 
distributions compared to the results from numerical simu¬ 
lations will help us further understand the range of physical 
conditions in galaxy halos. 


5 CONCLUSION 

We have investigated the equivalent width ratio 
lUMgii/WFeii for a large sample of Mgll absorption 
line syste ms from the JHU-SDS S Metal Absorption Line 
Catalog of IZhu &: MenardI (l2013l l and discovered that the 
ratio shows a trend with redshift. This trend does not 
appear to be the result of data artifacts or selection effects 
in the catalog. Several reasons may exist for the observed 
evolution. However, a simple model shows that the trend 
is naturally explained by the evolving Fe/Mg abundance 
ratio that results from the changing relative rates of 
Type la and core-collapse SNe. Star-formation (and hence 
supernovae production) is primarily located in galactic 
disks, whereas the absorption lines under investigation 
here form in circumgalactic clouds in galaxy halos at large 
galactocentric distances. Hence, if the equivalent width 
ratio reflects abundance variations in the CGM clouds, and 
if the enrichment by supernovae is indeed the primary cause 
for the CGM metallicity variations, then outflows must be 
efficient in transporting metals out of star-forming disks to 
large galactocentric radii. 

Future spectroscopic surveys with the next generation 
of highly multiplexed multi-object spectrographs (e.g., DESI 
on the Mayall telescopfQ, PFS on Subarru, and 4MOST on 
VISTAQ) will result in larger samples of absorption line sys¬ 
tems and the ability to explore the chemical evolution of 
the Universe in greater detail. Higher signal-to-noise ratio 
and higher resolution spectroscopy of large samples span¬ 
ning a large redshift range can confirm whether the trend 
observed in the Fe/Mg ratio is observed in other a-element 
lines and explore whether the changes are due to the chem¬ 
ical enrichment history of the Universe (as suggested here) 
or if evolution in the cloud properties (e.g., ionization, kine¬ 
matics, Doppler parameter, substructure, etc.) is instead 
responsible for the trend. Detailed comparisons between 
the data and state-of-the -art cosmological simulations (e.g., 
IVoeelsberger et 'al]|2014l l will be critical to this effort. 
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